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ABSTRACT 


*The  hypervlrlal  theorem  is  used  to  develop  expectation  value 
relationships  for  one  and  two  electron  central  force  problesis. 

Ground  state  helium  trial  functions  are  scaled  to  identically 
satisfy  a  particular  family  of  these  relationships.  It  is  shown  that 
this  procedure  results  in  improved  expectation  values  for  positive 


powers  of  the  radial  coordinate. 
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A  STUDY  OF  ONE  AND  TWO  El£CTRON  CENTRAL  FORCE  PROBLEMS 
USING  THE  HYPERVIRIAL  THEOREM 


Sidney  L.  Gordon 


The  hypervirial  theorein"*^  states  that 

(r,  c  H,wjr)  =  o 


(1) 


Here^  H  Is  the  Hamiltonian  operator. 

H'  is  a  bound  state  eigenfunction  of  H  • 

W  Is  called  a  hypervirial  operator,  and  Is  a  time  Independent 
function  of  the  coordinates  and  conjugate  momenta.  Vs/  must  be  chosen 
so  that  H  Is  hermit  ian  with  respect  to  'f'  and  W'f'  ,  e.  g. 

(v,  H  wr)  =  C  ,  (2) 

The  square  brackets  represent  the  conmutator  of  H  and  Vy  , 

defined  by  CH,Wj|  ^HNaZ-WH, 
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Epstein  and  Epstein  used  (1)  to  develop  expectation  value  relation¬ 
ships  for  the  non-relatlvlstlc  and  relativistic  hydrogen  atc'  ,  and  for 
the  Isotropic  harmonic  oscillator.  We  extend  their  procedure  to  the 
non-relatlvlstlc  one  and  two  electron  central  force  problems  in  sections 
I  and  II  of  this  paper,  and  develop  a  variety  of  hypervirial  relations. 

In  section  III,  we  scale  approximate  ground  state  helium  wave  functions 
to  satisfy  a  particular  family  of  these  relations. 


I.  One  Electron  Central  Force  Problem 

The  Hamiltonian  for  a  one  electron  central  force  problem  Is 
P*-  k** 

H  =  ^  Vcr) 

iLffi  zmr^ 


(3) 


where  t  Is  the  radial  coordinate  of  the  electron. 

p  is  the  momentum  conjugate  to  K"  ,  and  Is  explicitly  given  by 
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3  ® 

the  expression  p  •=  -  t  il  -Jr  ^  ( r  )  . 

Is  the  square  of  the  orbital  angular  momentum^  and  for  this 
problem  is  a  number  equal  to 

V  is  the  central  force  potential  energy  and  is  a  function  of  T 

only. 

Let  us  consider  as  the  hypervirlal  operator 


w  =  -f  p  w 

where  'f  Is  a  function  of  r  .  Using  the  commutator  relations 

CAB^CJ  Ai:B>Cj-h  CA>C3B  (5) 


and 


C  p.-fl  = 


(6) 


we  find  that 


<CH.4P1>  =  -i±<|fp‘>.^‘<^,p) 


=.  O 


(7) 


where  the  angular  brackets  stand  for  the  expectation  value  using 
a  bound  state  eigenfunction.  The  expression  (7)  is  sero  because  of  the 
hypervirlal  theorem.  Now^ 

<[:H,fil>  =  o.  (8: 

Expanding  equation  (8)  we  find  that 


(9) 
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Multiplying  eq.  (3)  on  the  left  by  and  taking  the  expectation  value, 

s>r 

we  deduce  that 


J.  <  r-*  2±\ 

amSar*^  /  =\sr/  zS  ^  5r  / 

-  <V  21  >  . 

^  sr  ' 

Substituting  (9)  and. (10)  into  (7),  and  rearranging  terms  we  find  that 

E<l^>- -<<"•*>) 


(10) 


(11) 


/We  now  specialize  (11)  by  letting 


4  =  r" 


(12) 


where  t1  Is  an  Integer.  We  show  in  the  appendix  that  condition  (2)  is 
satisfied  when 


n  >  I  -  z  Jl 


(13) 


where  I  is  the  angular  momentum  quantum  number  of  the  state.  Then, 
after  rearranging. 


nE<r"-'>+  -  ^*]  /  r  ■»-») 

=  n<Vr"->  +  ^ 


(14) 


Equation  (14)  shows  how  the  expectation  values  of  the  powers  of  T  are 
influenced  by  the  type  of  central  potential.  If  we  let  V  *  *” 
and  V  “  ^  we  obtain  Epstein's  relations^  for  the 

nonrelativistlc  hydrogen  atom  and  Isotropic  harmonic  oscillator, 
respectively. 
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H.  Two  Electron  Central  Force  Problem 


The  Hamiltonian  for  the  two  electron  central  force  problem  is 
given  by 

H  -  II  -h  li.  -I-  JLl.  ^  !S:2l  . -hV,  i-Vx+ViA 

zm  zm  zror,**  amrj* 


zrn 


(15) 


where 


Is  the  radial  coordinate  of  electron  1. 
is  the  radial  coordinate  of  electron  2. 
pl  Is  the  momentum  conjugate  to  l|  . 

|l^  Is  the  momentum  conjugate  to  . 


with  a  similar  expression  for  «Sl 
Is  a  function  of  I',  . 
\/f^  Is  a  function  of  . 
V,^ia  a  function 
and  2.  Because  of  theV 


of  ,  the  distance  between  electrons  1 

term^  and  k.|^  are  not  nuad>ers. 


As  the  hypervlrlal  operator,  consider 


W  =  ^p, 


<i«) 


where  ^  is  a  function  of  and  .  Corresponding  to  equation 

(7) ,  we  obtain 

<  C>.,W3>  = 


(17) 


It  Is  not  possible  to  express  this  relation  entirely  In  terms  of  the 
radial  coordinates.  For  example,  corresponding  to  eq.  (9)  we  find  that 

/»±Lp  \  =  ii/?!±+  ?21_  \ 

9r»^'/  z  \  9r,9  9r,9V»  / 

-  <  Sr.'^). 


(18) 


However^  If  we  let  f  be  a  f-jactlon  of  only^  equation  (17)  almpll- 

fies  to 


(19) 


We  now  specialize  (19)  by  letting 

f  =  r," 


(20) 


where  II  Is  an  Integer.  It  is  shown  In  the  appendix  that  condition 
(2)  is  satisfied  for  an  S  state  If 


nil,  (21) 


Equation  (19)  then  becomes 

JOl  <  r  ^  ^  o-3S 

zm  ^  '  Si;; -  ^  / 

+  n"'®)  =  i  <  (fr.‘  ''i+v»))> .  (22) 


We  shall  study  eq.  (22)  in  the  next  section. 

The  potential  energy  terms  for  the  helium  problem  are  given  by 


V, 

n  > 

Vi 

and 

=  "• 

2 

“5r  » 

(23) 

= 

Then^  when  n  Is 

equal  to 

one^  equation  (22)  becomes 
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t 


<  T,>  =  T  <  n  (  Ir.  I V.  +  V.x))> 
=  -  i  <  V,  +  i  V.x  > 


(24) 


where  T|  Is  the  kinetic  energy  of  electron  1.  Equation  (24)  la  the 
one  electron  vlrlal  theorem. 

S  state  helium  functions  may  be  completely  expressed  In  terms  of 


three  variables ^  such  as 


and  is  the  angle 


between  and 

these  coordinates  Is 

.a 


).  The  S  state  helium  Hamiltonian  in  terms  of 


H  ^  II  + 


) 

(25) 


+  V,  Vj,  •4*  V,x 


The  relations  already  derived  in  this  section  are  valid  for  this 
Hamiltonian  if  we  let 


(26) 


If  we  consider  coordinate  systems  such  aa  the  Itylleraas  Ti  ,  , 

Via  ''B  must  be  careful.  The  Hamiltonian  may  no  longer  have  a 

simple  dependence  on  the  radial  momenta.  In  r,  >  *Ti  >  coord- 

InateS)  the  partial  derivatives  with  respect  to  Ti  and  are 

subject  to  the  triangle  equality,  and  the  variations  are  no  longer  along 
the  radial  direction. 


III.  Scaling  Helium  Trial  Functions 

In  sections  I  and  II,  we  used  the  hypervirial  theorem  to  derive 
expectation  value  relations  for  one  and  two  electron  central  force 
hamiltonians.  Ife  now  change  our  point  of  view,  and  consider 

F  (  i)  »  (5,CH,w]5) 


(27) 
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as  a  functional  of  $  .  Then,  because  of  the  hypervlrlal  theorem,  if 
^  is  equal  to  an  exact  bound  state  eigenfunction  of  H  >  F  ( 
will  be  Identically  equal  to  zero.  If  ^  Is  equal  to  an  approximate 
trial  function  with  an  adjustable  parameter  K  >  e*g. 


i  =  <l>  ( K)  , 


(28) 


we  uay  determine  a  value  of  K  ,  say  Ko  >  such  that 

F  (  ♦(  Ko))  •=  o  .  (29) 

We  then  say  that  satisfies  the  hypervlrlal  theorem  correspond¬ 

ing  to  Vy  .  In  this  section  we  choose  K  to  be  a  scaling  parameter, 
and  adjust  approximate  ground  state  helium  functions  to  satisfy  the 
family  of  hypervirlal  theorems  generated  by 

W«  —  P»  n*  (30) 


The  helium  atom  hamlltonian  in  atomic  units  is 


^  Ik  2r,»  ar»*  r,  fx  Hx 


(31) 


The  corresponding  to  (30)  are  given  by  (22),  and  in  atomic 

units  are  (within  a  multiplicative  constant) 


+  2($.  (3 

-  i($,'""(i/v,  +  v.a.))$)  nx, 


He  choose  as  our  approximate  trial  function 

-Ky,„(rr,*rx) 


e 


K»,nCrrx+ 


(33) 
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where  normalization  constant; 

V  Is  a  radial  correlation  parameter; 

K  V'lV)  Is  the  scale  parameter,  and  is  determined  by  the 
condition  that 


Pn  (  4^y,n  )  ^  O  . 

Eq.  (34)  may  be  explicitly  solved  for  K  }f,n  giving 


(34) 


^rin  — 


(35) 


where  the  subscript  IT  on  the  angular  brackets  means  that  the  expec* 
tatlon  values  are  evaluated  using  the  unsealed  function 


-(vn  +  fi)  -(vnL+n) 

0  *t*  W 


ffe  therefore  generate  a  family  of  trial  functions  ^  If correspond- 
Ing  to  the  family  of  hypervlrial  operators  Wn. 

mien  V)  «  I  ,  our  procedure  corresponds  to  scaling  a  function 
to  satisfy  the  vlrlal  theorem,  and  therefore  (PlTtl  will  give  the 
lowest  energy  for  a  given  value  of  Y  mien  Y^  |  ,  ^  Y>  I 

is  equal  to  the  E chart -Kellner  function  *  ;  when  V  S  2.18/1.19  , 

^  y,  I  is  equal  to  the  Eckart-Hylleraas  radially  correlated  functlon^'^. 

In  figure  1,  we  have  plotted  K  YiO  versus  H  for  Y  ^  I 
and  Y  s  2. 18/1. 19.  The  values  for  T)  equals  1  thru  5  are  given 
In  table  I.  For  both  values  of  Y  ^  K  Y«V>  goes  asymptotically  to 
1  as  n  becomes  large.  The  curve  for  Y  s  2.18/1.19  ,  however, 
goes  thru  a  maxlsnim. 

A  series  of  curves  for  V  I  is  shown  In  figure  2.  For  all 
values  of  Y  ^  I  .  Ky,n  goes  asymptotically  to  1  as  V) 
becomes  large.  The  behavior  of  K  ViV)  when  0<  Y  I  isay  be 
deduced  by  noting  that 


^  >  n 


1 


K  V,  n 


V 


(36) 
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We  therefore  conclude  that 


Lifn 

n  -^oo 


K 


¥»T1) 


9  when  Y  1  I  ) 

X  when  o  <  Y  i  I , 
V 


(37) 


Physically,  this  means  that  for  large  n  and  this  choice  of  trial 
function,  the  hypervirial  tbecrems  generated  by  w„*  r,"p,  are 

asymptotically  satisfied  uhen  the  most  extended  orbital  becomes  a  hydro¬ 
gen  atom  orbital.  The  ccnvergence  of  (37)  is  more  rapid  the  more 
localized  the  least  extended  orbital. 

As  stated  above,  gives  the  best  value  of  the  energy  for 

a  given  value  of  V  .  However,  this  does  not  imply  that  4^  V) I  gives 

the  best  expectation  values  for  other  operators.  As  an  example,  we 

consider  the  positive  powers  of  the  radial  coordinate  of  an  electron  in 

helium,  r,^  . 

We  have  listed  in  table  II  values  of  ^  Ti  ^  for  ^  ~  1  thru  5  , 

which  we  have  calculated  using  the  Eckart-Hylleraas  function^'  ,  and  the 
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Eckart-Hylleraas  function  with  (9.p)  configuration  interaction  .  We 
have  also  given  the  accurate  values  for  C  and  ^  calculated 

by  Fekeris^^.  By  comparing  these  numbers,  it  seems  reasonable  to  assume 
that  the  true  values  of  <r.*>  are  somewhat  lower  then  the  values 
calculated  using  the  second  function. 

The  values  of  ^  calculated  using  Y  ■  1  are 

given  in  table  III.  Numbers  corresponding  to  ^  and  n  ranging  from 
1  thru  5  are  given.  The  results  for  Y  <■  2.18/1.19  are  given  in  table 
IV.  The  last  column  in  both  tables  III  and  IV  contains  our  estimates  of 
the  true  values  of  <r.*>  from  table  II.  The  first  columns  of 
expectation  values  in  tables  II  and  IV  are  identical  because 
is  equal  to  the  Eckart-Hylleraas  function  when  ■■  1  and  Y  >  2.18/1.19. 

By  comparing  these  numbers,  we  see  that  for  both  values  of  Y  , 
the  functions  corresponding  to  VI  »  2  give  the  best  overall  values  of 
<n-*>  .  This  is  true,  even  though  the  functions  corresponding  to 
n  ■>  1  give  the  best  value  of  the  energy.  The  results  are  particularly 
striking  for  V  «  1  .  The  Eckart-Hylleraas  function  gives  fairly 
reliable  results  for  the  first  5  positive  powers  of  ,  however  the 
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n  "  2  function  does  give  a  noticeable  Improvement. 

The  expectation  values  of  positive  powers  of  the  radial  coordinate 
depend  strongly  on  the  behavior  of  the  wave  function  In  the  region  of 
configuration  space  where  one  electron  Is  far  from  the  nucleus.  Tfe 
conclude,  therefore,  that  scaling  approximate  trial  functions  to  satisfy 
the  hypervlrlal  theorem  for  vv  *  P|  improves  their  quality  In 
this  region  of  configuration  space. 

In  contrast,  scaling  trial  functions  to  satisfy  the  hypervlrlal 
theorem  for  *  ITi  P|  Improves  their  quality  in  the  region  of 

configuration  space  where  both  electrons  are  close  to  the  nucleus.  We 
may  adjust  Vin  to  simultaneously  satisfy  the  hypervlrlal  theorems 
for  n  ■  1  and  2  by  treating  V  as  a  variable.  Then,  V  Is  deter¬ 
mined  by  the  condition  that 

Kt»i  =  Ky»a  .  (38) 

« 

We  are.  In  effect,  finding  a  value  of  V  which  allows  us  to  scale  ♦r.n 
to  simultaneously  satisfy  the  hypervlrlal  theorems  generated  by  W  x  Ti  p| 

and  r,**  p,  . 

Equation  (38)  Is  satisfied  when  V  ■  1.822  and  ■  1.195  .  The 
resulting  values  of  K,  are  given  In  table  V.  For  comparison,  we 

list  the  values  of  ^  calculated  using  when  Y  ■  2.18/1.19, 

and  n  ■  1  and  2  .  We  also  list  our  estimates  of  the  true  values  of 
from  table  11.  There  Is  only  a  small  difference  between  the 
values  of  ^  Tg^^  because  Y  -  2.18/1.19  ■  1.832  alsiost  satisfies 
(38).  The  calculated  values  of  ^  t  however,  are  best  when  trial 

function  (33)  Is  scaled  to  simultaneously  satisfy  the  first  and  second 
hypervlrlal  theorems. 

We  may  extend  this  procedure,  and  adjust  trial  functions  to  satisfy 
the  first  •  hypervlrlal  theorems  by  including  •  adjustable  param¬ 
eters.  It  Is  reasonable  to  expect  further  improvements  In  the 
expectation  values  of  <  Tg^  .  However,  to  verify  this,  we  need  reliable 
values  for  expectation  values  of  higher  powers  of  the  radial  coordinate 
In  helium.  We  are  presently  learning  how  to  obtain  these  ntimbers,  and  ^ 

will  report  the  results  of  these  Investigations  at  a  later  date. 


11 


Acknowledgment 

The  author  acknowledges  with  thanks  the  many  helpful  discussions 
with  Professor  J.  0.  Hlrschfelder  during  the  course  of  this  work. 


Appendix 


12 


•f 


In  this  appendix  we  derive  the  limitations  on  ITI  given  by  the 
Inequalities  (13)  and  (21).  Acceptable  hypervlrlal  operators  must 
satisfy  equation  (2).  Let  us  define 

■0=  Cf,  HWV)  -  (  H'I'.Wf )  .  (39) 

Equation  (2)  then  becomes 


!>  =  O.  (^0) 

We  use  atomic  units  In  the  following  discussion. 

One  Electron  Central  Force  Problem 
The  hamlltonlan  Is  given  by 

H  ^  VCr),  (41) 

Substituting  (41)  into  (39)  we  find  that 


Singularities  may  occur  at  T  *  0  ;  therefore^  we  enclose  the  origin 
with  a  spherical  surface  S  of  infinitesimally  small  radius  C  ,  and 
define  X)  by  Its  principal  value.  Therefore, 

T>  =  -T  Lim  f  (43) 

*  •'♦o  It' 

where  means  that  the  Integration  is  performed  over  all  space  except 
for  the  region  enclosed  by  the  spherical  surface  S  .  We  may  now  apply 
Green’s  Theorem  in  the  Second  Form,  then 


A  llm  ^  r*»'* 

^  t4o  1$ 


(44) 
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The  sign  has  changed  because  we  cunsider  the  normal  to  the  surface  to 
be  pointing  away  .•rom  the  origin. 

Equation  (44)  reduces  to 

T)  =  i  Lim  €*f  (45) 

^  t  *♦©  I  r  •"  c 

where  the  Integration  is  ever  the  solid  angle.  However 

-  K<r)  '©-(©, 4>)  ,  (46) 


and  therefore 


We  let 

W  =  <-"0  =  4-r*-'i.fr  )  .  (48) 

I  dr 

Substituting  (48)  into  (47)  and  expanding  we  obtain 


V  - 


\  L;m  r  <«-.)  R*-*"  n  R  «"♦' 

+  (  r«!!:l  -  (i&f 

A  Njr/  Jr*€* 


(49) 


We  assume  that 


r^Vtr)  =Q  , 

r o 

Then^  near  the  orlgin^^ 

a  =  Ar^ 


(50) 


(51) 


where  A  is  a  constant,  and  oCT^) stands  for  an  expression  which 
satisfies  the  condition  that 
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L  I W  O  C  -  C  >  <52) 

r-^o  r* 

where  -C  is  a  constant  (which  may  be  zero).  Substituting  (51)  into 
(49)  we  find  that 


T>  = 


Zi 


e-*o  *“ 


(53) 


Therefore,  D  ■  0  when  n  ,  orV>'il"3L-X  .  This  proves 

inequality  (13). 


W 


Two  Electron  Central  Force  Problem  (  S  state) 

The  hamiltonian  is  given  by 

H  =  V,cr.)+  VxCro  v»a.cr,a) 


(54) 


For  an 

S  state 

,  is  real. 

D  = 

S| 

T  Li'w 

We  let 

Vs/ 

Then,  i 

sorresponding  to  (49) 

D  = 

-D,  +  TJi 

(55) 


(56) 


(57) 


where 
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T5.  =  X 
2  1 

Lim 

fe  -♦o 

-f  ( 

and 

Lim 
e  -»o 

i)  +  n  >1'  83: 

»*■. 


For  an  3  state^  if 


L  im  r,^  V,  <''■•)  =  O  , 

r,  -^o  (58) 


L  im  Vx  C^x)  =  O  , 

fx-^O 

and 


Li^  ^  constant  which  may  be  sero^ 

r,x->o 

then  V  is  a  non  zero  finite  fimction  of  Tx  when  Tf  is  zero,  and 


vice  versa. 

LVj  X 

If  we  also  assume  this 

to 

be  true  of  the  partial 

derivatives 

of 

V  ,  then 

D.  = 

J_ 

0 

21 

f 

fr-*©  ^ 

A  <»»-») 

-8* 

and 

1 

• 

2( 

Li  m  $ 
8  ->o  1 

a'  C  + 

o 

<*"♦»)] 

^  • 

Here^  A  and  are  non  zero  constants. 

Therefore^  "D  -  0  for  an  5  state  if  FI  ^  |  ^  which  proves 

inequality  (21). 
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Table  L 


n  as  calculated  from  equation  (35). 


n 

l.n 

^2.18/1. 19 

1 

1.688 

1.191 

2 

1.500 

1.193 

3 

1.344 

1.266 

4 

1.227 

1.444 

5 

1.145 

1.744 

Table  II. 

^  y  as  calculated  from  three  helium  ground 
state  wave  functions. 


<'i'>  ‘ 

n  atomic  units 

t 

Bckart-Hylleraas^'^ 

Eckart-Hylleraas 
with  (2p)2  C.I.’ 

»  u  j  10 

Pekerls 

1 

0.94 

0.93 

0.9295 

2 

1.24 

1.22 

I. 1935 

3 

2.14 

2.10 

4 

4.66 

4.53 

5 

12.28 

11.81 

Table  III 


as  calculated  from  with  "V  »  1  .  Each  column  gives 

values  of  <'i'>  corresponding  to  a  definite  value  of  n  .  The  value 
of  n  is  given  at  the  top  of  the  column.  The  last  column  gives  our 
estimates  of  the  true  values  of  <'i'>  from  table  II. 


in  atomic  units 

n 

t 

1 

2 

3 

4 

5 

estimated 

1 

0.89 

1.00 

1.12 

1.22 

1.31 

0.93 

2 

1.05 

1.33 

1.66 

1.99 

2.29 

1.19 

3 

1.56 

2.22 

3.09 

4.06 

5.00 

<2.1 

4 

2.77 

4.44 

6. 90 

9.94 

13.11 

<  4.5 

5 

5.76 

10.37 

17.97 

28.37 

40.10 

<11.8 

Table  IV. 

as  calculated  from  with  T  2.18/1.19  ■  1.832  .  The 

arrangement  of  this  table  is  the  same  as  table  III. 


<'i'> 

in  atomic  units 

n 

t 

1 

2 

3 

4 

5 

estimated 

1 

0.938 

0.936 

0.88 

0.77 

0.64 

0.930 

2 

1.24 

1.23 

1.09 

0.84 

0.58 

1.19 

3 

2.14 

2.13 

1.78 

1.20 

0.68 

<2.1 

4 

4.66 

4.63 

3.64 

2.16 

1.01 

<  4.5 

5 

12.28 

12.17 

9.03 

4.68 

1.82 

<11.8 

L9 


Table  V. 

<'i'>  as  calculated  from  V^f)  with  V*  ■■  2.18/1.19  and  1.822. 
The  last  column  gives  our  estimates  of  the  true  values  of  from 

table  II. 


In  atomic  units 


y  -  2.18/1.19  -  1.832 


0.938 

1.236 

2.143 

4.66 

12.28 


y  -  1.822 


n  »  2 

n  ■  1  and  2 

estimated 

0.936 

0.936 

0.930 

1.231 

1.231 

1.19 

2.130 

2.128 

<2.1 

4.63 

4.61 

<  4.5 

12.17 

12.11 

<11.8 

